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Abstract 

We employ the recently developed worldline numerics, which combines string-inspired field 
■ theory methods with Monte Carlo techniques, to develop an algorithm for the computation 

of pair-production rates in scalar QED for inhomogeneous background fields. We test the 
algorithm with the classic Sauter potential, for which we compute the local production rate for 
the first time. Furthermore, we study the production rate for a superposition of a constant E 
field and a spatially oscillating field for various oscillation frequencies. Our results reveal that 
the approximation by a local derivative expansion already fails for frequencies small compared 
' to the electron mass scale, whereas for strongly oscillating fields a derivative expansion for 

, the averaged field represents an acceptable approximation. The worldline picture makes the 

f"^ ■ nonlocal nature of pair production transparent and facilitates a profound understanding of 

l/^ ' this important quantum phenomenon. 

o ' 
in 

O ; 1 Introduction 

O |' Pair production was first proposed for electron-positron pairs in strong, temporally and spatially 

■ , constant electric fields El E] ■ Today it is often referred to as the Schwinger ^ mechanism. 

As a nonperturbative mechanism, pair production is of great theoretical interest. From a phe- 
nomenological point of view, it corresponds to probing the theory in the domain of strong fields. 
Consequently, we encounter pair production in many topics of contemporary physics, for instance, 
black hole evaporation [5] and e + e~ creation in the vicinity of charged black holes |SJ[7] as well 
as particle production in hadronic collisions [8] and in the early universe EH ■ Since QED pair 
production in strong fields represents the conceptually simplest case, it can serve as a theoretical 
laboratory for all these cases. 

A sizeable rate for spontaneous pair production requires extraordinary strong electric fields, 
comparable in size to the so-called critical field strength, which corresponds to the electron-mass 
scale, E CT = m 2 /e f=a 1.3 • 10 18 ^-. For a long time, it seemed inconceivable to produce macroscopic 
electric fields of the required strength in the laboratory, but today, with the development of 
strong lasers, there are several promising experiments in progress ^JEIEJ; for a discussion of 
experimental requirements, see |14| . 

Many different theoretical methods, such as the propertime method WKB techniques 

1171 1181 IT§| , the Schrodinger- Functional approach [2U| , functional techniques [221 > kinetic 
equations [23 El 123 ESI , various instanton techniques E3 EH EH E0| , Borcl summation [SUES, 
and propagator constructions |3HI l3"4"| , have been developed to study pair production in external 
fields. Also, finite-temperature contributions have been determined which first occur at the two- 
loop level jSSl EE] ■ Of particular conceptual interest is the production rate in terms of the effective 
action for a given background, which is also used in this work. Owing to an intimate relation 
between the effective action and the vacuum-persistence amplitude, it is the imaginary part of 
the effective action that encodes information about pair production which, in this context, is 
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interpreted as spontaneous vacuum decay. This approach yields the instantaneous production 
rate, neglecting back-reactions and memory effects. However, this rate can serve as a source 
term for kinetic equations, which can then take back-reactions and memory effects into account 

[231 EH [23 urn 

Even though the existing methods follow a well defined and technically stringent concept, their 
application often faces serious technical and conceptual difficulties. Up to now, no reliable and 
universal method be it analytic or numeric — is available for the calculation of pair-production 
rates in inhomogeneous electric fields. In standard approaches, functional traces have to be eval- 
uated with the knowledge of the spectrum of the corresponding differential operator, which is 
only available for special cases. Moreover, controlling the divergencies that possibly occur upon 
summing up the eigenvalues is a delicate task. 

In the present work, we solve these problems by using the recently developed numerical world- 
line techniques j33 QEl l-ffll 1401 l4~T] which are based on the string-inspired worldline formalism 
[4*21 |4~*"1 14*11 l4l"l IT71 14*51 """""1 l*"T"] . The important advantage compared to other approaches lies in 
the fact that worldline numerics can be formulated independently of any symmetry of the back- 
ground. The identification of and the summation over the spectrum of quantum fluctuations are 
done in one single and finite step. For simplicity, we confine ourselves to scalar QED; generalization 
to spinor QED is, in principle, straightforward and will be discussed below. 

Beyond the computational advantages of worldline techniques, the worldline picture also helps 
to understand conceptual aspects in more depth. In particular, the nature and the role of nonlo- 
calities become highly transparent from the worldline viewpoint, since the worldlines themselves 
represent extended virtual trajectories of the fluctuating particles in coordinate space. In the 
present context, we are aiming at the quantum effective action which, of course, receives nonlocal 
contributions in general. However, many standard approximation methods suppress (or shade) 
nonlocalities by construction, as, e.g., the derivative expansion. Hence, pair production as de- 
scribed by the Schwingcr formula is often recognized as a nonperturbative phenomenon, but not 
so much as a nonlocal phenomenon. Nevertheless, the latter property is crucial, as the following 
heuristic argument elucidates: in order for a virtual pair to become real, i.e., on-shcll, the pair 
must gain at least the amount of 2m of energy; this is only possible by propagating in opposite 
directions in the electric field. This derealization of the pair wave function is mandatory for 
gaining sufficient energy. 

In constant electric fields, this derealization remains invisible in the final result. By contrast, 
in inhomogeneous fields the spacetime dependence of the delocalized wave function matters a 
great deal and can even dominate the resulting effect, as our results demonstrate. In the worldline 
picture, the nonlocal effects already become transparent on the level of the formalism, since the 
extended worldlines exactly describe the derealization of a virtual pair. 

At this point, we would like to stress the difference of the present work to earlier applications 
of worldline numerics. Whereas the algorithms developed so far in 1881 IT-iT?! 1401 141j have proven 
their capabilities for computing the real part of the effective action (and action densities), the 
computation of the imaginary part is by no means a straightforward generalization. The reason 
for this lies in the truly Minkowskian nature of the problem of pair production: vacuum decay 
only occurs for real, i.e., Minkowskian, electric fields. This contrasts with the indispensable 
necessity of a Euclidean formulation for solving the worldline integrals by a statistical Monte 
Carlo algorithm. In practice, this results in an overlap problem: the finite Euclidean worldline 
ensemble can have little overlap with those worldlines that contribute dominantly to Minkowski- 
valued observables. We solve this fundamental problem by resorting to a technique developed in 
jo 1 " in the different context of nonperturbative Euclidean worldline numerics: we fit a suitable 
cumulative density function (CDF) of the Euclidean ensemble to a physically motivated ansatz that 
can be continued analytically to Minkowski space. We should emphasize that this continuation 
represents an extrapolation of certain ensemble properties to Minkowski space which is an a 
priori uncontrolled procedure resulting in systematic errors. We check this extrapolation carefully 
against various analytically known results and find negligibly small systematic errors compared to 
the statistical Monte Carlo errors. Hence, we regard the overlap problem as solved for the present 
problem. This solution is obtained at the expense of numerical cost; moreover, the algorithm 
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can, in principle, not be made arbitrarily precise, in contrast to former applications of worldline 
numerics. Nevertheless, for the problem of pair production and as far as the experimentally 
required accuracy is concerned, we believe that our algorithm is sufficiently powerful. 



2 Worldline formalism for pair production 

The vacuum-persistence amplitude can be related to the effective action Tm in Minkowski space, 

(n\e" lHT \Q) = e lFM . 

The corresponding probability for the vacuum to decay spontaneously is 

p = \ — g-2imr M 

In the case of QED with electric background fields, vacuum decay occurs in the form of spontaneous 
pair production, the production rate per unit time and volume of which is directly proportional 
to the imaginary part of the effective action density (effective Lagrangian). 

In scalar QED, the one- loop contribution to the Euclidean effective action Te reads 

r E L4] =lndet(-(9 + ieA) 2 + m 2 ), (1) 

where Tm and Te differ by a minus sign, Tm = — Te- In the worldline approach, the logarithm of 
the determinant in _D-dimensional spacetime is represented by a path integral |50j . 

where the integration parameter T is called the propcrtimc. The path integral runs over all closed 
worldlincs, parameterized by the propcrtimc. The worldlines can be viewed as the trajectories of 
the virtual fluctuations in coordinate space. The path integral is normalized to give 1 in the limit 
of zero gauge potential. We split the path integral into an integral over all paths with a common 
center of mass xq and an ordinary integral over all xq, x(t) — * xo + x(t), where dr x(t) = 0. 
Introducing the Wilson loop, 

W X0 [ X {t)\ := er le So dT iA{x +x{r)) ; (3) 

and its expectation value, 

(W X0 }:=[ Vx{T)W X0 [x(r)]e-^ dT ^, (4) 



(0) = x(T) 
CM 



we can write: 



^ - -J^yJJl J d°-o I ^j~ 2 e- 2T (W X0 ) + c.t. (5) 

Here wc have added counterterms (c.t.) which have to be fixed by renormalization of physical 
parameters. If the electric field is nonzero, r 1 [A] obtains an imaginary part, arising from poles of 
the Wilson- loop expectation value (W Xo ) on the real T axis. Surrounding the poles by halfcirclcs 
in the upper half plain in agreement with causality leads to 

ImTUA] = -—^uH I d ° Xo Im E rr, ! +D/2 ^ m2T ^(-m)Res((W X0 ),T pol ) , (6) 

( 4?r ) J T pol 7 P°1 

where the sum goes over all poles with positions T po i- The exponential factor with Gaufiian velocity 
weight in the path integral in Eq. suppresses the contribution of long paths. Therefore, the 
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integral is dominated by paths that tightly wiggle around the common center of mass. This 
gives rise to the picture of a loop cloud sitting at xq and scanning the background field in the 
neighborhood of xq. Hence, the nonlocal nature of the phenomenon is already apparent in the 
formalism. 

Let us mention in passing that the path integral for a constant E background is Gaufiian, can 
thus be done exactly, and results in (W Xo ) = eET/ s'm(eET); see below. Summing over the pole 
positions of the inverse sine results in the famous Schwingcr formula (for scalar QED in this case), 

ImT^E = const.] = -_ {eEf £ i-J- e"* ™ (7) 

n— 1 

displaying the nonperturbative dependence on eE 0|; here, V denotes the space-time volume. 
Each term in the sum corresponds to production of n coherent pairs. 



3 Worldline numerics 
3.1 Worldline discretization 

The worldline numerical algorithm for the present problem partly resembles closely those developed 
in detail in jS3 IH81 1891 1401 RT] , the essential steps of which we will recall in the following for 
completeness. As a first step, we introduce the unit loop y(t), 

y(t) := -J=*(Tt). (8) 



The Wilson-loop expectation value then reads 

(W X0 ):= f r Vy(t)W Xo e-fo<*4. (9) 



»(0) = »(1) 
CM 



The exponential velocity distribution is now independent of T, whereas the Wilson loop yields 

W x = e~ ie f° dt ^Ty A {^ T y+ x °) . (10) 

For a finite ensemble of paths that are distributed according to the weight factor cxp(— J* dr \), 
the Wilson-loop expectation value is equal to the arithmetic ensemble average of W xo . Since the 
weight factor for unit loops is independent of T and xq , such a loop ensemble has to be generated 
only once for computing (W Xo ) for different T and xq. 

For the numerics, we discretize the propertime such that each loop y(t) is represented by a 
finite number of N points per loop (ppl) yk at t = k/N, with k = 1, • • • , iV. 1 We use the vloop 
algorithm 0U] to create an ensemble of m, discrete and closed unit loops {yk} with the distribution 
functional 

^[{2/fe}] = '5(yi + --- + yJv)exp^-^^(y fc - 2 / fc _ 1 ) 2 ^ , (11) 

with the condition yo = y^ for closed loops. Eq. 111|) represents the discrete form of the weight 
factor cxp(— J Q dr ^-), with the delta function reflecting the center-of-mass condition. 

For a gauge-invariant discretization of the Wilson loop, the gauge field, in principle, should be 
treated as a link variable, i.e., dty(t) A(y(t)) — > (yk — Vk-i) A((yk + yk+i)/2), with the gauge field 
evaluated at the center of the links. However, the link centers do not carry the same information 
about the distribution of the worldlincs in spacetime as the sites yk do: for instance, the link 
centers have a smaller average distance to the center of mass xo than the sites do. The use of the 

1 The worldline points yj. live in continuous space-time, S R D . For an alternative lattice formulation, see 
l52ll53l . 
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link centers actually corresponds to effectively shrinking the loop cloud. Of course, this difference 
becomes irrelevant in the propertime continuum limit TV — * oo. However, for small TV, this effect 
leads to sizeable systematic deviations from the continuum limit. We avoid this systematics by 
evaluating the gauge field at the sites instead, 

/ dtyA(Vfy + x )^y2(yk+i -y k )A(VTy k +x ). (12) 
J ° k=i 

It turns out that possible violations of gauge invariance for smooth gauges such as the Lorcnz 
gauge remain much smaller than other systematic and statistical errors for the background fields 
studied in this work. 2 

For the effective action and the pair-production rate, the T integration in Eq. nas to be 
performed. For the simple case of a constant field, this can be done elegantly by a fast-Fourier 
transform (FFT) after the T integration has been rotated onto the imaginary axis. Thereby, the 
pair production is obtained for a whole spectrum of masses and field strengths, respectively, all 
at once. This procedure and its limitations will be discussed in the Appendix. However, for 
more general field configurations, an overlap problem arises: when performing the T integration, 
one faces situations in which the path integral is dominated by very elongated loops, despite 
the exponential suppression by the weight factor. Physically, those virtual pairs that delocalize 
strongly gain more energy and have a larger probability of becoming real. In this case, the 
finite loop ensemble with only a few elongated worldlincs is no longer representative for the over- 
countably many paths of the path integral. To solve this problem, we have developed the routine 
presented in the following. Its cornerstone is a probability distribution analysis of particular 
worldlinc-cnscmblc properties along the lines suggested in |51) . 



3.2 CDF fit for pair production 

In order to motivate our algorithm, let us first consider the case of a constant homogeneous 
electric field E in Minkowski space. 3 This is related to the Euclidean gauge potential by A\& = 
(0,0,0, — iExi) T . The corresponding Wilson loop can be written as 



W{I) = e- TeEI , where I := [ 

Jo 



dty^yx. (13) 



The scalar quantity / contains all relevant information about the unit loop for the present case. 
The probability density function (PDF) of / for our loop ensembles is defined by 



(j-J dty A y^j e-ti dt €. (14) 



P(I)= / VyS\I 

J CM 

With the aid of a Fourier representation of the 6 function, the path integral becomes Gaufiian and 
yields 

P(I) = J cosh- 2 (|/) (15) 
for constant fields. In terms of the PDF, the Wilson-loop expectation value can be written as 

/oo 
dl P(I)W(I), (16) 
-oo 

resulting in (W) = TeE / sm{TeE) in agreement with the Schwinger pair-production rate for 
constant fields, cf. Eq J7J). 

2 In the general case, we, of course, recommend the gauge-invariant link variable discretization. In order to 
reduce the systematic error mentioned above, order 1/JV improvements of the action may be useful. 
3 E denotes the Minkowskian, i.e., physical, field strength. 
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For inhomogeneous field configurations, (W Xo ) can be computed in a similar way 
the definition of I, 

_ i /p dt yA(VTy + x ) 
1x0 - jTEo ' 

the PDF becomes space-time and proper-time dependent, 

Px (I)= I Vy8{I-I X0 )e-ti dt %. (18) 

/y(0)=»(l) 
CM 

But for each space-time point, the Wilson-loop average can still be computed analogously to 
Eq. (HI, with W(I) = e~ TeEo1 similar to Eq. JT^. The reference field strength E is a priori 
arbitrary and has been introduced to obtain a dimcnsionless quantity. In most cases, we may 
use the local field strength Eq := \E(xo)\, or some averaged value. For the constant E field, 
our generalized definition of I Xo conforms to the previous one. The PDF of I Xo is generally not 
known analytically but will be computed numerically from a finite loop ensemble. Nevertheless, 
analytical knowledge about P Xo (I) is required, owing to the following reasons: 

• The use of a Monte Carlo algorithm does not only demand the worldline spacetime metric 
to be Euclidean, but also requires the contour of the propertime integral to run along the 
real T axis. However, as is already obvious for the constant-field case, the integral in Eq. 

is well defined only for \TeE\ < ir. At \TeE\ = tt, the first pole T pol of (W Xo ) is hit. For 
larger values of \eET\, the I integral has to be replaced by its analytic continuation, which 
can only be constructed if P Xo (I) is known analytically. 

• By using finite loop ensembles, we already face an overlap problem for small T values: the 
majority of loops have a small I value, whereas those few loops with large I dominate 
the I integral in Eq. (|16|l : see the Appendix. A controlled extrapolation of the PDF to 
large / values from reasonably big worldline ensembles can thus reduce the numerical cost 
considerably. This can be achieved by fitting the numerical PDF data to an analytical 
ansatz. 



Generalizing 
(17) 



The last point is, of course, related to the nonlocal features of pair production. The quantity / 
on the one hand is connected to the electrostatic energy gain of a virtual pair that propagates in 
a background field, and on the other hand roughly measures the space-time extent of a worldline. 
The dominance of large / values in the final result arises from strongly dclocalized virtual pairs. 

To obtain an analytical estimate for the PDF, we generalize the result for the constant field, 
Eq. (|15|) . by the following ansatz governed by two parameters a and v: 

P XB (I) = N cosh" 21 ' go/). (19) 

The parameters control the two main features of the distribution: width and sheerness. Both 
parameters depend on the spacetime point xq and on the propertime parameter T. The normal- 
ization constant iV is a function of a and v fixed by J dl P Xo (I) = 1. Numerically more convenient 
is the corresponding cumulative density function (CDF) of |7|, 



r\I\ 

D X0 (\I\)= / dIP(I) 
J-\i\ 



For given values T and xq, we determine a and v by a fit of the numerical data to this CDF. 
Inserting the resulting parameters into Eq. I|19(l yields the desired analytical expression for P Xo (I). 
Performing the / integral in Eq. (|16|l gives (W Xo ) as function of a and v, 

W ) — N - 7ra na_J_ 

{Wxo) - ira T{2v) 
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This result also represents the desired analytical continuation to arbitrary values of T or | e^o^l ■> 
and solves the problem of Wick rotating the result of the Euclidean path integral back to Minkowski 
space. We observe that the second Gamma function in the numerator is responsible for the pole 
structure of (W Xo ) on the positive real axis. Poles occur if 

T V P 

v — = -Z, with/ = 0,1, 2,-- . (20) 

7ra 

Since a and v depend on T, Eq. i|20|) determines the pole positions T po i only implicitly; in practice, 
we solve for T po \ iteratively. At the pole location, the corresponding residue is 

v ' dT \ rra ) 

which we plug into Eq. © to obtain the pair-production rate. For a reliable control of the 
statistical error, we perform a jackknife analysis for all secondary quantities. For the systematic 
error due to the propertime discretization, we increase the number of N ppl to approach the 
continuum limit at least within the statistical errors. 

Obviously, the reliability of our results depends crucially on the ansatz (|19|l for the PDF. Apart 
from our consistency arguments referring to the shape of the PDF and the resulting pole structure, 
final support can only be given by nontrivial tests described below. In summary, the sufficiency 
of the ansatz is confirmed by the following arguments: 

• The free parameters control the two most essential features of the distribution, the width and 
the sheerness, which encode, in particular, the important contributions from the strongly 
delocalized virtual pairs. 

• The exact functional form for the constant-field limit is supported by the ansatz. Even 
without further checks, we could thus expect satisfactory results at least for slowly varying 
fields. 

• As a nontrivial analytical confirmation, we stress that the ansatz leads to a reasonable pole 
structure of (W) that can encode information about coherent n pair production. 

• The ansatz provides highly convincing results for the Sauter potential including the constant 
field limit as special case, as presented in the next section. Any systematic deviations from 
the exact result are negligibly small compared to the statistical error. 



4 Sauter potential 

The Sauter potential defines an electric field with solitonic profile in one spatial direction which 
is constant in all other directions including time. The direction of the field vector is constant and 
coincides with the solitonic-profile direction. An analytical expression of the corresponding total 
pair-production rate has been found by Nikishov |54j . In Minkowski space, the Sauter potential 
reads 

A a \ M = -atanh^x 1 ), A l \ M = 0, E 1 ^ = ° 

cosh (kx 1 ) 

The parameter k defines the inverse width of the electric field, whereas a governs its maximum, 
£ mal[ = ak. The constant-field limit is recovered for k — > for fixed ak. 

As an example, Fig. ^ shows the x\ dependence of the local pair-production rate ~ Im£ ff 
for k = 0.4m and E max = ak = (m 2 /e) computed by our algorithm. It is compared to the 
approximated effective Lagrangian obtained by a derivative expansion to lowest order, i.e., by 
assuming the field to be locally constant (Schwingcr formula). We observe that the local rate 
predicted by the algorithm is spatially smeared compared to the Schwinger formula. 
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Figure 1: Spatial distribution of the effective Lagrangian's imaginary part for a Sauter potential. 
The numerical result is compared to the locally constant-field approximation that overestimates 
the true result by up to ~ 50%. Parameters of the Sauter potential: k = 0.4m, E max = (m 2 /e). 
Parameters of the loop cloud: n L = 100000, N = 1000 ppl. 
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Figure 2: The imaginary part of the effective Lagrangian in the center of a Sauter potential with 
maximal field strength m 2 /e versus the inverse width parameter k. The dashed line marks the 
analytically obtained contribution of the first pole to |InxC e ff| f° r the constant-field limit (higher 
poles give a 1% correction). n L = 100000, N = 1000 ppl. 



The pair-production density in the center x 1 = of the Sauter potential with maximal field 
strength E mSj ^ = (m 2 /e) is shown in Fig. Inversus the width parameter k; units are set by the 
electron mass scale. For large width k — > 0, the constant-limit is approached, and our CDF 
fit algorithm correctly reproduces the Schwinger formula. The more interesting limit occurs for 
k = m where the production rate vanishes. Even though the electric field is still nonzero, the 
width of the Sauter potential is equal to the Compton wavelength. Therefore, even if a virtual 
pair delocalizes completely along the direction of field lines with the e~ going to a; 1 — ► oo and the 
e + going to — oo, the pair cannot acquire enough energy to become real. This important physical 
example is missed completely by the locally constant-field approximation, emphasizing the role of 
nonlocalities. 

Moreover, the limiting case of k — > m is an extreme and crucial test for our algorithm based 
on the PDF ansatz Ijl9(l : in the vicinity of this limit, there is literally not a single worldline 
in our finite ensemble that exhibits the strong derealization required for giving rise to a direct 
contribution to the final result (the number of sufficiently elongated worldlines is exponentially 
suppressed). Nevertheless, the overall distribution of / values allows for a controlled extrapolation 
via the CDF fit, leading to a numerical estimate even for the directly inaccessible regime. As a 
measure for the resulting error, we mention that our result for the case k = m is not exactly zero, 
but |Im£ e ff|/m 4 = 5.73 • 10 -8 ± 1.03 • 10 -6 . We conclude that possible systematic errors induced 
by our CDF fit algorithm are negligibly small compared to the statistical error. 

Finally, Fig. shows the integrated total pair-production rate Imr compared to the Nikishov 
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result. The agreement is satisfactory and the vanishing pair production for ea = to is reproduced 
within the error bars. 
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Figure 3: The imaginary part of the effective action for a Sauter potential as fraction of the locally 
constant-field approximation Imri c versus the width parameter k in units of to: comparison of the 
numerical result with Nikishov's analytic expression, = 100000, N = 1000 ppl. 



5 Sine-modulated potential 

In this section, we study the superposition of a spatially varying sine potential with a constant field. 
This configuration is of general interest, as it is representative for a class of field configurations 
which are superpositions of a slowly varying field — in our example the constant field — and higher- 
oscillation modes. A very important aspect is the dependence of the pair-production rate on the 
spatial oscillation frequency of the small-scale field structures. We consider this example as a 
paradigm for the role of nonlocal phenomena in pair production. 
In Minkowski space, the potential is given by 

A°\ M = -asinfW) - Eqx 1 , A i \ M = 0. 

It corresponds to an E field in x 1 direction with field strength 

£^|m = E a + afccos(fcx 1 ), 

which has extremal field strength of -Emax. mm = Eq ± ak. As an example, we study a field with 
E Q = 0.2(TO 2 /e) and £ raax = 0.3(m 2 /e). 

Figure 0] shows the position of the first pole T po \ of the Wilson- loop expectation value on the 
real propertime axis for x 1 in the center of a maximum of the field strength. For small fc, the 
pole position of the constant-field limit E = E max is reproduced. For large k, the pole position 
converges to the result of the averaged field E = Eq. In between, the curve is not monotonically 
increasing, as one might have expected, but reaches T values which are significantly larger than 
in both limiting cases. As a consequence, the corresponding local production rate will be smaller 
than in the constant-field limit E = Eq. 

This behavior is, of course, a consequence of nonlocalities and can be easily understood in the 
worldlinc picture in terms of loop clouds: Starting with the limit k — > 0, a loop cloud sitting at a 
maximum detects a constant field of strength E max . A sketch of this scenario if given in Fig. EJa). 
If k is increased and the wavelength of the sine becomes shorter, the loop cloud overlaps more and 
more with the minima on cither side of the maximum and the pole moves to larger T values. If k 
exceeds a certain value, in our example at about k = 0.8m, the two minima close by dominate the 
Wilson-loop expectation value, Fig. Efb). Despite the maximum in the center of the loop cloud, 
the pole is at a larger T value than for the averaged field. Not until the loop cloud approaches 
the adjacent maxima, Fig. |5fc). do the T values become smaller again, finally converging to the 
value of the averaged field, Fig. Eld). 
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Figure 4: Position of the first pole T po \ of (W) on the real proper-time axis at a maximum of the 
field strength. With increasing frequency k, the pole moves from the constant-field limit E = E mayi 
to the limit E = Eq. Parameters of the field: Eq = 0.2(m 2 /e), E max = 0.3(m 2 /e). In between, 
it develops an unexpected maximum corresponding to a minimum of the local production rate. 
Parameters of the loop cloud: n L = 100000, N — 1000 ppl. 




Figure 5: An artist's view on a loop cloud (worldline ensemble) at a maximum of the field strength. 
For small frequencies, it detects only the maximum (a). After increasing the frequency, the two 
nearest minima dominate (b). For larger frequencies the cloud encounters further maxima (c), 
until it perceives an averaged field (d) . 

Since the Wilson-loop expectation value at a maximum of the field strength can be dominated 
by the adjacent minima, the inverse situation can also occur at a minimum where the result can be 
dominated by the two adjacent maxima. In this case, the first pole of (W) is at a smaller T value 
than for the averaged field, leading to a larger imaginary part of the effective Lagrangian. This 
inversion is shown in Fig. where the spatial distribution of the imaginary part of the effective 
action for k = 1.8m is plotted in comparison to the constant-field limit E = Eq. We observe that 
the nonlocalitics induce a seemingly paradoxical phenomenon in this case: the maxima of the local 
pair-production rate occur at the minima of the electric field strength and vice versa. 

Figured depicts the imaginary part of the total effective action per space-time volume for our 
example configuration versus the frequency k. In contrast to its density at Xo, ImT does not fall 
below the result for the averaged field. For oscillation frequencies near k = 0, we observe that 
the locally constant-field approximation based on the derivative expansion fails rather early by 
an order of magnitude for k ~ 0.5m; this is remarkable, since the effective expansion parameter 
k 2 /m 2 ~ 0.25 might have been considered as small enough. 

In the opposite limit, for large frequencies k, we obtain the averaged constant-field limit E = 
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Figure 6: Spatial distribution of the imaginary part of the effective-action density for the sine- 
modulated potential with k = 1.8m compared to the constant-field limit E = E$. Nonlocal effects 
lead to the seemingly paradoxical phenomenon that the pair-production rate is maximal at the 
field-strength minima and vice versa, ni, = 200000, TV = 1000 ppl. 
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Figure 7: The imaginary part of the total effective action per space-time volume against the 
frequency k. The dashed lines mark the locally constant-field approximation and the result for 
the averaged field E = Eq, respectively. The former (dashed lines) misses the true result by an 
order of magnitude already for k/m ~ 0.5. til = 200000, N = 1000 ppl. 



Eq. It is remarkable that the imaginary part of the effective action reaches the value of the 
averaged field for k values as small as about k = m, whereas its density still fluctuates spatially 
for even larger k values, as seen in Fig. The fluctuations cancel each other, so that they have 
no effect on the integrated quantity. The numerical accuracy does not eliminate the possibility of 
a A;-dependent structure for k values larger than m. According to the values of Fig. |7| the central 
values suggest a slight increase of the pair production for k > m, until it falls back to the result 
for the averaged field if k/m — » oo. To definitely clarify this, larger loop ensembles are necessary 
at the expense of CPU time. However, the present result shows that any possible k dependence 
for k > m has to be relatively small and the averagcd-field approximation yields good results in 
this range. 

Let us finally compare our results for the spatially sine-modulated field with those for spatially 
homogeneous fields with time dependencies. Especially the case of an electric field oscillating in 
time with frequency w has been studied with WKB methods ^3 El 1181 Ej which were originally 
developed for ionization processes in atomic physics |55| . The nature of pair production in this case 
depends on the size of the "adiabaticity parameter" 7 := muj/(eE); for small 7 <C 1, the result 
approaches the Schwinger formula and pair production thus is a nonperturbative phenomenon. 
For large 7 > 1, the result becomes perturbative in (eE) / (mu>) and pair production arises from 
multiphoton scattering. In our case, we can, of course, also form a similar parameter 4 7 = mk/eEo, 
with 7 small or large roughly corresponding to the two limiting cases discussed above. However, it 

4 This parameter is not unique in our case, since the parameter a of the sine modulation introduces yet another 
scale. 
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is important to stress that pair production is nonpcrturbative in both limits for our sine-modulated 
field. In particular, the l&rge-k/m (or large 7) case cannot be understood in terms of multiphoton 
processes. Taking the external field to all orders into account is essential for the final result. 

6 Conclusion 

We have developed a new universal approach for computing local production rates for spontaneous 
pair creation by the Schwingcr mechanism in scalar QED. Our method is based on the combination 
of the worldlinc formalism with Monte Carlo techniques. As a first result, we have not only 
rediscovered Nikishov's analytic result for the total pair-production rate in a Sauter potential, but 
moreover we have computed the local pair-production rate for this classic case for the first time. 
Most importantly, the algorithm is not restricted to any spatial symmetry of the given background 
potential but is applicable for arbitrary potentials. 

As a nontrivial example, we have applied the algorithm to a constant electric field modulated by 
a spatial sine oscillation. This field configuration is representative for a whole class of fields with 
large-scale structures and small-scale oscillations. By varying the spatial oscillation frequency, 
qualitatively different features of pair production can be investigated. For small frequencies, our 
numerical result agrees with the derivative expansion to lowest order; the latter breaks down 
completely for spatial variations on the order of a few times the Compton wavelength. On this 
length scale and below, our results show clearly that another approximation scheme becomes 
reliable: the local production rate can well be approximated by inserting the spatially averaged 
field into the Schwingcr formula. This averaged-ficld approximation can be trusted on the few- 
percent level for spatial variations of the size of the Compton wavelength. We would like to 
emphasize that the small validity bound of the derivative expansion for the imaginary part of 
the effective action density is not related to the same observation for the real part, as discovered 
in [321; the latter arises from a subtle interplay between nonlocal quantum contributions and 
local counterterms, whereas the imaginary part is not affected by renormalization counterterms. 
Furthermore, the derivative expansion for the real part of the integrated effective action works well 
even for Compton-scale variations |5fij . whereas it breaks down early for the imaginary part, as 
displayed in Fig. [7| 

Apart from these quantitative results for the particular field configurations considered here, our 
findings emphasize the crucial role of nonlocalities in the phenomenon of pair production. Without 
the feature of delocalization of a virtual pair, spontaneous vacuum decay would not occur. The 
worldlinc picture underlying our algorithm is particularly powerful in capturing these nonlocalities 
and, moreover, understanding their consequences in an intuitive way. Especially our results for 
local pair-production rates illustrate the nature and the role of nonlocalities transparently. For 
instance, the seemingly paradoxical situation that maxima of pair-production rates can occur 
at minima of the field strength (cf. Sect. EJ cannot be understood from a local approximation. 
However, the worldline picture identifies a natural explanation of this phenomenon in terms of the 
delocalization properties of the virtual pairs described by the worldline trajectory. 

From a technical perspective, we have developed a numerical Monte Carlo algorithm that on 
the one hand requires a Euclidean formulation for the quantum fluctuations, but on the other hand 
produces reliable results for truly Minkowski-valucd physical observables. The inherent overlap 
problem is solved in the present context by a physically motivated ansatz for a suitable cumulative 
distribution function (CDF) to which the numerical data can be fitted and that can be analytically 
continued to Minkowski space. Even though the success of this procedure depends strongly on 
the problem at hand, we believe that such techniques can be useful in other Minkowski-valucd 
problems as well. The algorithmic strategy itself has been invented in the context of Euclidean 
field theory [51j . where it has turned out to be highly powerful in a study of nonperturbative 
worldline dynamics. 

Several extensions of our work are desirable and possible. So far, we have only considered 
spatial inhomogencities, but any realistic field configuration will also exhibit variations in time. 
In fact, timelike variations bring in a new complication, since our Monte Carlo worldlines live in 
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imaginary time, whereas physical fields depend on real time. Therefore, our algorithm is directly 
applicable to all those cases where the physical field is known analytically, such that its analytic 
continuation to imaginary time can be evaluated and plugged into the numerics. For instance, the 
exact result for a solitonic profile in time direction as solved in [37] will be a benchmark test for 
such an investigation. 

Furthermore, our results can, in principle, straightforwardly be generalized to ordinary spinor 
QED. As a new complication, the Pauli term ~ cr^ v F^ v occurs in the worldline integrand. Since 
this term depends also on the worldline trajectory, the probability distribution function (PDF) 
of the ensemble will not only depend on the quantity I as defined in Eq. 1|17[) . but also on the 
worldline averaged Pauli-term exponential; let us denote the latter with J, which is also a scalar. 
Our algorithm might be generalized as follows: first, compute the PDF of J from the ensemble 
and bin the loops according to their J value. Then, apply the present algorithm to each J bin 
separately; in particular, the same analytic-continuation technique can be used. Finally, integrate 
over J with the aid of the PDF of J. It is important to note that the J integral can be done last, 
since the Pauli-term worldline average cannot induce any poles for the proper-time integral. Of 
course, since each relevant J bin has to contain sufficiently many worldlincs, this generalization of 
our algorithm will at least be an order of magnitude more time consuming than the one for scalar 
QED. At this point, we should stress that the computations for the present work have still been 
performed on ordinary desktop PC's. 

Finally, it is instructive to compare our method to the instanton technique of [271 I3U| , where 
the instanton approximation of the worldline integral has been shown to give the leading-order 
contribution to pair production. For instance, in the constant-field case, the one-pair-production 
rate is generated by one instanton which is a circular loop. Small fluctuations around this path 
lead to the correct imaginary prefactor. In comparison to this, our worldlines are extraordinarily 
complex. Not a single worldline loop in our ensembles resembles a circle or fluctuations thereof. 
This gives rise to the conjecture that the computation of the imaginary part requires very little 
information about the shape of the loops. We expect that we should be able to extract the 
instantonic content of our loops by a suitable cooling procedure that removes large-amplitude 
fluctuations. In view of the success of the instanton approximation, only instantonic plus small- 
amplitude-fluctuation information appears to be relevant for pair production. This agrees with our 
observation that pair production is induced by delocalized "large" loops that can acquire enough 
energy in the E field. Therefore, it is well possible that a different loop discretization which 
optimizes instantonic properties allows for an even more efficient computation of the imaginary 
part. A further investigation of this topic may lead to an even deeper understanding of pair 
production. 
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A Constant field: straightforward approach 

In the following, we present a straightforward realization of worldline numerics for calculating the 
pair-production rate in the constant-field case. The algorithm presented here is an immediate 
generalization of the standard algorithm successfully used for the real part of the effective action 
[S3 EH1 CM EH GI] For a constant field in four dimensions, Eq. © reads 
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where I is defined as in Eq. (|13J) . and the counterterms for on-shell renormalization are included. 
Rotating the T integration contour onto the imaginary axis and substituting s = —iTeE yields a 
Fourier integral, 

r ^(f) 2 r$-*/^(<-"'>- i+ r)- 

If the worldlinc-cnsemble average (e~ lIs ) can be computed reliably, Eq. I)22J1 offers a highly efficient 
algorithm with the aid of the FFT: in this case, F^. can be computed for a whole spectrum of 
frequencies m 2 /eE all at once with FFT. The resulting imaginary part is shown in Fig. |H1 It 
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Figure 8: Imaginary part of the effective action obtained by FFT. til = 5000, N = 1000 ppl. 

is highly remarkable that this numerical procedure gives satisfactory results in a wide range of 
scales, extending over five orders of magnitude, with little consumption of CPU time. However, 
the algorithm fails for small field strengths. The precise limit is given by the size of the largest 
loops in the finite loop ensemble: only loops with |7| values larger than m 2 /eE contribute to the 
imaginary part of T^. For weak fields, this implies that only a few or even no loops contribute 
and the computation fails. 

Beside this problem which is already relevant for the constant-field case, there is a second 
limitation. For a different contour in the complex T plane which supports large ReT values, the 
Wilson-loop expectation value is dominated by the loop with the largest I value. The Monte Carlo 
algorithms break down here, since the error bars become as large as the central value. In general, 
for inhomogeneous background fields, it is not possible to find a suitable integration contour to 
avoid this problem. 

These limitations of the straightforward approach are a manifestation of the fact that the 
Euclidean worldline ensemble has insufficient overlap with Minkowski- valued obscrvablcs for weak 
fields. 
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